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Abstract 

A multi-linear variable separation approach is developed to solve a differential-difference 
Toda equation. The semi-discrete form of the continuous universal formula is found for 
a suitable potential of the differential-difference Toda system. Abundant semi-discrete 
localized coherent structures of the potential can be found by appropriately selecting the 
arbitrary functions of the semi-discrete form of the universal formula. 
PACS numbers:02.30.Ik, 05.45.Yv 

1 Introduction 

It is very difficult to find explicit exact solutions both for nonlinear partial differential equations 
(PDEs) and for nonlinear differential-difference equations (DDEs). In linear mathematical physics, 
the Fourier transform and the variable separation approach are two of the most effective ways to find 
exact solutions of linear equations. In nonlinear mathematical physics, the so-called inverse scattering 
transform (1ST) has served as "nonlinear" Fourier transform for nonlinear integrable models. However, 
it is very hard to extend the variable separation approach to nonlinear cases even for integrable systems. 
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Recently, several kinds of "variable separation" approaches, say, the classical method, the differen- 
tial Stackel matrix approach [1], the geometric method [2], the ansatz-based method [3, 2], functional 
variable separation approach[4], the derivative dependent functional variable separation approach[5], 
the formal variable separation approach (nonlincarization of the Lax pairs or symmetry constraints) 
[6] and the multilinear variable separation approach (MVSA) [7]-[9] are in progress. Among these 
approaches, the MVSA may be the most powerful method to solve many (2+l)-dimensional systems 
both for integrable models[7]-[9] and for nonintegrable systems[10]. In [9], a quite universal formula 

jj _ 2(aia 2 - apa^gyPx ^ 
_ (a + aip + a 2 q + a 3 pq) 2 ' 

is established to describe suitable physical quantities for various (2+l)-dimensional models solvable 

via MVSA. In (1), a , a\, a 2 and a 3 are arbitrary constants and p is an arbitrary function of {x, t} 

for all of the known MVSA solvable models, while q of (1) may be an arbitrary function of {y, t} 

for some of MVSA solvable models, or an arbitrary solution of a Riccati equation for some other 

MVSA solvable models. Because some arbitrary characteristics, lower dimensional functions (like p), 

have been included in the universal formula (1), by selecting them appropriately, abundant localized 

structures like the multiple solitoffs, dromions, lumps, breathers, instantons, ring solitons and chaotic 

and fractal patterns have been found[9]. 

In this paper, we are interested in the following important question: Can the MVSA be extended to 

solve some nonlinear DDEs? In section 2, the MVSA is outlined for arbitrary DDEs. In sections 3 and 

4, the MVSA is applied to a special differential-difference Toda equation (SDDTE). A semi-discrete 

form of the universal formula (1) for a suitable potential of the SDDTE is given in section 5. Starting 

from the semi-discrete form of the universal formula, the abundant semi-discrete localized excitations 

can be found for the potential of the SDDTE. Some special semi-discrete localized excitations of the 

model are also plotted in section 5. Section 6 contains a short summary and discussions. 

2 Outline of the MVSA 

Consider a given DDE of the following type 

F(x i ,n J ,u(x i ,n j - k, ...),u Xi (xi,n,j - k, ...), = 1, ...,Ni,j = 1, ...,N 2 ,k = 0,±1, ±2, ...) 
= F(u) = 0, (2) 

where Xi, i = 1, 2, ...Ni are continuous variables and rij, j = 1, 2, ...N2 are discrete variables. As 
in the continuous case, the MVSA can be performed for DDEs via the following standard procedures. 
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(i) . Multilinearize the original DDE (2) by using a suitable Backhand transformation. Usually, the 
resulting equations are the bilinear equations for integrable systems. 

(ii) . Choose a seed solution of the Backlund transformation as general as possible with one or more 
arbitrary functions. 

(iii) . Make a suitable variable separation ansatz with some variable separated functions. Usually, the 
ansatz is just the generalization of the Hirota's two-soliton solution. 

(iv) . Substitute the ansatz into the multilinear equations and separate the resulting equations to 
several variable separated ones. 

(v) . Solve the variable separated equations. Usually, to solve the variable separated PDEs and/or 
DDEs is still very difficult for any fixed seed solution. However, one can treat the problem in an 
alternative easy way: The variable separated functions appeared in the ansatz can be considered as 
arbitrary functions and then the function(s) appeared in the seed solution should be fixed from the 
variable separated equations. 

(vi) . Find a suitable field or potential which possesses a special variable separated solution described 
by the universal formula (1) for PDEs (7V 2 = in (2)) or a semi-discrete form of (1) for DDEs. 

(vii) . Discuss the possible semi-discrete localized excitations by selecting the arbitrary functions 
appropriately. 

To see the details on the procedures of the MVSA for the DDE systems, we take the SDDTE as 
a simple concrete example in the remained sections. 



In nonlinear discrete and semi-discrete physics systems, the most famous and important systems are 
the so-called Toda systems which are widely used in physics [11]. In this paper, we only consider a 
special differential-difference Toda equation (SDDTE) 



Q(n) yt = exp[Q(n + 1) - Q(n)][Q(n + 1) + Q(n)] y - exp[Q(n) - Q(n - l)][Q(n) + Q(n - 1)]„, (3) 



where Q(n) = Q(n, y, t) is a function of the discrete variable n and the continuous variables {y, t}. 
The SDDTE (3) was firstly derived by Cao, Geng and Wu in a remarkable paper [12]. Some interesting 
integrable properties of the SDDTE (3) have been given in [12] and [13]. 

To solve the SDDTE (3) via MVSA, one can use the following dependent variable transformation 



3 SDDTE and its generalized bilinear form 




(4) 
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to bilinearize it. In the transformation (4), p(n) = p(n,t), the seed solution of the SDDTE has been 
selected as an arbitrary function of {n, t} for convenience later. 

Substitution of the dependent variable transformation (4) into the SDDTE yields 



, 7i _ , } { DyDtf(n) ■ f(n) - 2exp[p(n) - p(n - l)]D v exp(D ra )/(n) ■ fjn) ^ = () 



2/(n)2 ' " - (5) 

where T + is a shift operator, i.e.,T + F(n) = F(n + 1), and the Hirota's bilinear differential operator 
D™D% and the bilinear difference operator expD„ are defined by 

k 



D™D k t a ■ b 



(- 


_o_y 








dy'J 




dt'J 



cxp(SD n )a(n) ■ b(n) = exp 



a(n)b(n) 



v=y',t=t' 
= a(n + 5)b{n — 5). 



d__J)_ 

s dn dn' 

Multiplying (5) by the inverse operator of T + — 1 leads to a generalized bilinear SDDTE 



D y D t f(n) ■ f(n) - J(n, t)D y exp(£»„)/(n) • f(n) + R(y, t)f(n) 2 = 0, 



(6) 



where 



J(n) = 2exp[p(n) — p(n — 1)] 



(7) 



is an arbitrary function of {n, t} and R(y,t), the kernel of the difference operator T + — 1, is an 
arbitrary function of {y, t}. 



4 Variable separation solution of the SDDTE 

In order to find some exact solution of (6), similar to the continuous cases[7]-[9], we look for the 
solutions of (6) in the form 

f(n) = a + aip(n, t) + a 2 q(y, t) + a 3 p(n, t)q(y, t), (8) 

where ao, &i, a 2 and as are arbitrary constants and the variable separation functions q(y,t) = q and 
p(n, t) = p(n) are only functions of {y, t} and {n, t} respectively. (8) looks like the Hirota's two 
soliton form when q and p(n) are exponential functions. Substituting the ansatz (8) into (6), we have 

-2(a 2 + a 3 p(n)) 2 q y qt + 2(a 3 a - a 2 a 1 )q y p nt + 2(a 2 + a 3 p(n))(a + aip(ra) + a 2 q + a 3 p(n))q yt 
-J(n)(p(n + 1) -p(n - l))(a 3 a - a 2 ai)q y + R(a + a 2 q + a\p(n) + a 3 qp(n)) 2 = (9) 
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Because p(n) and J(n) are only functions of {n, t} and q and R are only functions of {y, t}, (9) can 
be separated into the following two equations, 

qt = (ao + a 2 q) 2 ci + (ai + a 3 q) 2 c 2 + (a + a 2 q)(a 1 + a 3 q)c 3 , (10) 
p t (n) = (a a 3 - aid - 2)(c 2 - c 3 p(n) + cip(n) 2 ) + ^ J(n)(p(n + 1) - p(n - 1)), (11) 

when the arbitrary function R is selected as 

i? = -2{c\a% + c 2 aj + c 3 a 2 a 3 )q y , (12) 

with Ci = Ci(t), c 2 = c 2 (t) and c 3 = c 3 (t) being arbitrary functions of t. 

In principle, as long as the arbitrary functions a, c 2 , c 3 and p(n) (and then J(n)) are fixed, we 
can obtain the corresponding special solutions of the (10) and (11) and then the special solutions of 
the SDDTE (3). However, it is still very difficult to solve the nonlinear DDE (11) for fixed nonzero 
J(n). Fortunately, as in the continuous cases discussed in [7]-[9], because of the arbitrariness of 
the function J(n), we can treat the problem alternatively. We can consider the function p(n) as an 
arbitrary function of the variables n and t and fix the function J(n) from the equation (11). The 
result reads 

2 

J (n) = — — — [(a a 3 - a x a - 2) (c 2 - c 3 p(n) + ap(n) ) - p t (n) . (13) 

p{n - 1) - p(n + 1) 

It should be pointed out that the Riccati equation (10) is totally same as that of the asymmetric 
Nizhnik-Novikov-Veselov (ANNV) equation[9]. To find out some special solutions of (10), one may 
select the arbitrary functions appropriately. Here we list two special selections. 
(1). If we write c\, c 2 and c 3 as 

a\A 2t a 3 (ai + a 3 A 2 )A lt _ (ai + a 3 A 2 ) 2 A 3t 

(aia 2 -a a 3 ) 2 (a x a 2 - a a a 3 ) 2 A x {a x a 2 - a a 3 ) 2 Ai ' 

c _ a 2 A 2t a 2 (ao + a 2 A 2 )A lt _ jap + a 2 A 2 ) 2 A 3t 

(aia 2 -a a 3 ) 2 (aia 2 - a a 3 ) 2 Ai (aia 2 - a a 3 ) 2 Ai ' 



c 3 = 



(a a 3 + a x a 2 + 2a 2 a 3 A 2 )A lt 2a 2 a 3 A 2t 



(a\a 2 - a a 3 ) 2 Ai (aia 2 - a a 3 ) 2 

(ao + a 2 ^4 2 )(ai + a 3 A 2 )A 3t 
(aiai-AfAx 

with A\ = A\(t) 1 A 2 = A 2 (t) and A 3 = A 3 (t) being arbitrary functions of t, then the general 
solution of (10) with (14)-(16) reads 
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where F\ = F\(y) is an arbitrary function of y. 
(2). If we select c\, c 2 and c 3 as 

ajbot a 3 (ai + a 3 b )b lt [jai + «3^o) 2 - b\al]b 2 t 

(aia 2 - a a 3 ) 2 (aia 2 - a a 3 ) 2 6i (<na 2 - a a 3 ) 2 6i 

^ _ flj&ot a 2 (flo + a 2 b )b lt _ [(a + a 2 fe ) 2 - albj]b 2t ^ 

(aia 2 - a a 3 ) 2 (aia 2 - a a 3 ) 2 6i (<na 2 - a a 3 ) 2 6i 



c 3 



3 a 3 + aia 2 + 2a 2 a 3 b )b u 2a 2 a 3 h 



0/. 



(oia 2 - a a 3 ) 2 &i (aia 2 - a a 3 ) 2 

, [(ao + ft2^ o)(ai + « 3 M - Q 2 a 3 ^f]^ 2 t 
(oia 2 - A) 2 6x 



+ 2 TTT ^2! ( . ZU - ) 



with b = b (t), b\ = bi(t) and b 2 = b 2 (t) being arbitrary functions of t, then the general solution 
of (10) with (18)-(20) reads 

q = Manh(& 2 + F 2 (y)) + & (21) 
with F 2 = F 2 (y) being an arbitrary function of y. 



5 Abundant coherent structures for a potential of the SDDTE 

Substituting all the results of the last section into (4), one can get many kinds of exact solutions for 
the field Q of the SDDTE. In continuous cases, it has been pointed out that for every of the MVSA 
solvable systems listed in [9], there exits a quantity that can be described by the universal formula 
(1). Now an important question is: 

Is there a suitable potential for the SDDTE that can be described by a suitable semi- 
discrete form of the universal formula (1)? 

Fortunately, it is straightforward to prove that if we define a potential of the SDDTE as 

u = -2Q y (n), (22) 

then 

u = U(n) = 2g y (a 2 ai - a 3 a )(p(n + 1) - p(n)) 

(ao + a 2 q + aip(n) + a 3 qp(n))(ao + a 2 q + aip(n + 1) + a 3 qp(n + 1)) 

It is clear that the function U(n) defined in (23) is just one suitable semi-discrete form of the continuous 
universal quantity U given in (1). We say the function U(n) is semi-discrete means that it is discrete 
in one direction and continuous in other direction. 
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Now starting from the semi-discrete form of the universal quantity, we can obtain abundant semi- 
discrete localized excitations for the SDDTE by selecting the arbitrary functions appropriately. 

The detailed studies show us that the semi-discrete localized structures for the potential u ex- 
pressed by the semi-discrete form of the universal quantity are very similar to the continuous ones 
which have been discussed in Ref. [9]. So in this paper we will not discuss all the possible localized 
excitations but only list some particular examples. 

Example 1. Resonant semi-discrete dromion and solitoff solutions. 

If we restrict the functions p(n) and q of (1) as 

N N 

P( n ) = ^exp(fcin + ujit + x 0i ) = ^exp(^) (24) 

i=l i=l 
M J 

q = ^2 ex P(^y + Vot) ^2 exp(Oji + t j), (25) 
i=\ j=i 

where xoi, yoi, toj, ki, w,, Ki and tti are arbitrary constants and M, N and J are arbitrary positive 
integers, then we have a single resonant semi-discrete dromion solution or semi-discrete multiple solitoff 
solutions. The selection (25) is related to the selections on the functions of Ai, i = 1, 2, 3, Fx of (17) 

as 

A 3 =A 2 =0, (26) 
j 

A x =Y,eMttjt + t 0j ), (27) 
i=i 

Fl = ^i L T V (28) 

E,=i exp(^iy + yoi) 

and the c i; i = 1, 2, 3 are given by (14)-(16) with (26) and (27). 

In Fig.l, we plot four typical structures caused by the resonant effects of four straight-line semi- 
discrete soliton solutions. 

Fig. la shows the structure of a first type of single resonant semi-discrete dromion solution shown 
by (22) with (24), (25), 

M = N = 2, J = h = Kx = 1, k 2 = K 2 = i, oo = l, a 1 = a 2 = 10, a 3 = ^ (29) 

and 

xoi = Voi = hi = x 02 = J/02 = (30) 

at t = 0. 
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Fig. lb is a plot of a single resonant semi-discrete solitoff solution shown by (22) with (24), (25), 
(30) and 

M = N = 2, J = ki = K x = 1, k 2 = K 2 = \, a = 1, a ± = a 2 = 3, a 3 = 0, i = 0. (31) 

Fig. lc shows the structure of a second type of single semi-discrete resonant dromion solution 
shown by (22) with (24), (25), (30) and 

M = N = 2,J = k 1 = -K x = 1, -k 2 = K 2 = \ a = a 2 = 10, a 3 = 1, a x = \, t = 0. (32) 

Fig. Id is a plot of a four solitoff solution shown by (22) with (24), (25), (30) and 

M = N = 2, J = h = -Kx = 1, -k 2 =K 2 = \, a = 1, a x = a 2 = 3, a 3 = 0, t = 0. (33) 

Example 2. Semi-discrete oscillating dromions and lumps. 

If some periodic functions in space variables are included in the functions p(n) and q of (22), we 
may obtain some types of semi-discrete multi-dromion and multi-lump solutions with oscillating tails. 
The oscillated lump solution plotted in Fig. 2 is related to 

q = l + [y(cos(y) + 5/4)]2' p{n) = 1 + (n - ctf ' (34) 
«o = «3 = 1, o,i = a 2 = 5 (35) 

at t = 0. 

Example 3. Multiple ring soliton solutions. 

In high dimensions, in addition to the point-like localized coherent excitations, there may be some 
other types of physically significant localized excitation. Recently, we have found some different kinds 
of ring soliton solutions which are not identically equal to zero at some closed (2+l)-dimensional and 
(3+l)-dimensional curves and decay exponentially away from the curves[14, 9, 15]. 

In Fig. 3, a typical saddle type semi-discrete ring soliton solution is plotted for the potential u 
with the selections 

q = 6XP (~80 + 5 ) ' P{n) = 6XP ( {U 80 Ct) ) ' (36) 

and 

ao = «3 = 0, ai = a 2 = 5, (37) 

at t = 0. 
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Fig.la 




Figure 1: Four typical semi-discrete structures of SD0TE for the potential u expressed by (22) with (24) and 
(25). (a). A special single-peak semi-discrete resonant dromion solution, (b). A single semi-discrete solitoff 
solution, (d). A multi-peak semi-discrete dromion solution, (b). A semi-discrete four-solitoff solution. 



Fig.2 




Figure 2: Plot of a special oscillating lump solution of the SDDTE for the potential u expressed by (22) 
with (34) and (35) at t = 0. 



Fig.3 




Figure 3: Plot of a typical single saddle type semi-discrete ring soliton solution for the potential u of the 
SDDTE with the selections (36) and (37) at t = 0. 
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Because of the existence of the arbitrariness in the expression (22) for the potential u, there exist 
many kinds of multiple semi-discrete ring soliton solutions. Since the situations are quite similar to 
those in the continuous cases, we do not discuss them further. However it is worth to mention that 
as pointed out in Refs. [9] and [10], the interactions among the travelling ring soliton solutions are 
completely elastic without phase shift. 

6 Summary and discussion 

In the previous studies [7] -[9], we have successfully solved several famous (2+l)-dimensional non- 
linear continuous integrable models via a multilinear variable separation approach (MVS A). In this 
paper, we have further extended the MVSA to solve nonlinear differential-difference systems. Taking 
a special differential-difference Toda equation (SDDTA) as a concrete example, we have successfully 
solved the model via the MVSA. 

In continuous cases, a quite universal formula has been found to describe suitable fields or po- 
tentials of MVSA solvable models. For the SDDTA, a semi-discrete form of the universal formula is 
found for a suitable potential. An arbitrary function related to the discrete space variable is included 
in the semi-discrete form of the universal formula. Another function included in the formula is an 
arbitrary solution of the Riccati equation and the Riccati equation is totally same as that in some 
continuous MVSA solvable models. 

By selecting the arbitrary functions appropriately, one can find abundant semi-discrete localized 
excitations like the multiple solitoffs, dromions, lumps, breathers, instantons and ring soliton solutions. 
The semi-discrete localized solutions are quite similar to those of continuous cases shown in [9] . 

Though the MVSA have been successfully applied to the SDDTE, there are various interesting 
and important problems are worth to studying further. For instance, in continuous case, more than 
ten models have been solved by the MVSA (see Ref. [9] and the references therein) . How many DDEs 
may be MVSA solvable? It is also known that for a continuous integrable model, there exist some 
different types of integrable discrete forms. Is the semi-discrete form of the universal formula unique? 
In other words, how universal the semi-discrete quantity (23) is? In addition to the DDEs, there may 
be various fully discrete systems in real physics. Can we extend the MVSA to solve some fully discrete 
nonlinear systems? 

The authors are in debt to the helpful discussions with Drs. X-y Tang, C-l Chen and S-l Zhang. 
The work is supported by the Outstanding Youth Foundation (No. 19925522), the Research Fund for 
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